AD-A1Q3  270  ABERDEEN  UNIV  (SCOTLAND)  DEPT  OF  ENGINEERING  F/G  1/3 

APPLICATION  OF  NUMERICAL  METHODS  TO  THE  CALCULATION  OF  ELECTROS— ETC (U) 
AUG  81  J  R  SMITH.  P  LEES.  D  MCALLISTER  AFOSR-80-0223 


UNCLASSIFIED 


EOARD-TR-81-9 


Al03270 


Grant  Number  ATOSR-80-0223 

APPLICATION  OP  NUMERICAL  METHODS  TO  THE  CALCULATION  OF  ELECTROSTATIC 
FIELDS  IN  AIRCRAFT  FUEL  TANKS 


J.R.  Smith. 
Department  of 
Uriivers itv  of 


Engineering, 
Aberdeen . 


F.  Lees  D.  McAllister 

Department  of  Engineering  Dept,  of  Engineering 
University  of  Aberdeen.  University  of  Aberdeen 


REPORT  DOCUMENTATION  PACE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  Report  Number 


2.  Govt  Accession  No,j3.  Recipient  *  Catalog  Number 


4.  Title  (and  Subtitle)  / 

Application  of  numerical  methods  to  the 
calculation  of  electrostatic  fields  in 
aircraft  fuel  tanks. 


7.  Author  (s>^i  /smith/'  P.  Lees 
D.  McAllister 


8.  Contract  or  Grant  Number 
•/ AFOSR-80-0223 !  , 


10.  Program  Element,  Project,  Task 
Area  &  Work  Unit  -Numbers 


AFOSR 

61102F 

2301/DI 


9.  Performing  Organization  Name  Addrca-S- 

Department  of  Engineering,  |  ,)  {.  j  .  I 
Univeri sty  of  Aberdeen, 

Scotland 

/  ‘ 


11.  Controlling  Office  Name  ami  Address 

European  Office  of  Aerospace  Research  ( 
and  Development 
Box  14,  FI’O  New  York  09510 


14.  Monitoring  Agency  Name  and  Addveca 

European  Office  of  Aerospace  Research  ( 
and  Development 

Box  14,  FPO  New  York  09510  ( 


I  16.  4.  I/.  Distribution  Statement 

Approved  for  public  release;  distribution  unlimited. 


AFWAL 

62203F 

3048/05 


13.  Number  of  Faxes 
22 


Hi. 

Supplementary 

Notes 

19. 

Ki’.  y  o  I  d  s 

Electrostatic 

fields,  finite  element  method;  fuel  tanks 

70.  Abstract 

\x 

>The  solution  of  electrostatic  field  problems  occurring  during  the 
refuelling  of  aircraft  fuel  tanks  containing  explosion  suppressant 
foams  is  discussed.  A  computational  model  of  a  fuel  tank  is  set  up 
and  tun  finite  element,  method  is  used  to  calculate  the  electrostati 
potential  distribution  within  the  tank. 

Accession  F^r 
NTIS  FT  \iicl  ~~ 

*  DTIC  T ‘  8  □ 

Uimnae- □ 
Ju.'.  t  if  icc  t  ion _ 


mm;  i  h./j 


By_ _ 

M-yrtTitr:  typT/ - 

Aval lot: 11 tr  Ccd< 
■■■  il.  and/or 
.  pecial 


This  report  has  been  reviewed  by  the  EOARD  Information  Office  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At 
NT IS  it  will  be  releasable  to  the  general  public,  including  foreign 
nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


OWEN  MANCARELLA 

Lt  Colonel,  USAF 

Director,  Aeronautical  Systems 


^  "/f 


GORDON  L.  HERMANN 
Lt  Colonel,  USAF 
Deputy  Commander 


ELECTROSTATIC  HAZARDS  IN  AIRCRAFT  FUEL  TANKS  DURING  REFUELLING 

The  generation  of  high  levels  of  electrostatic  charge  in  hydro¬ 
carbon  fuels  (such  as  JP-4)  during  aircraft  refuelling  has  long  been 
recognized  as  an  explosion  hazard.  The  fuel  is  pumped  at  high  flow 
rates  through  pipes,  hoses,  and  filter/separators,  and  hence  is 
exposed  to  relatively  large  liquid/solid  interfaces.  The  double  layer 
created  at  these  interfaces  coupled  with  the  movement  of  fuel  across 
them  leads  to  a  net  unipolar  charge  being  acquired  by  the  fuel  as  it 
is  swept  along.  The  high  charging  tendency  of  a  fuel  such  as  JP-4, 
coupled  with  a  low  conductivity  (typically  <  10  pS/m)  can  lead  to 
hazardous  charge  accumulation  in  the  receiving  tank.  If  resulting 
local  electrostatic  fields  on  the  fuel  surface  exceed  the  breakdown 
value  for  the  vapour  space,  electrostatic  discharges  may  occur.  Such 
a  discharge  may  be  incendive  if  (i)  it  has  sufficient  energy,  and 
(ii)  the  fuel/air  mixture  lies  in  a  flammable  range  (i.e.  will  support 
combustion) .  It  is  thus  patently  of  importance  to  be  able  to  estimate 
electrostatic  potential  and  field  distributions  in  fuel  tanks  for 
given  charge  distributions. 

The  introduction  of  polyurethane  foam  into  fuel  tanks  to  act  as 
an  explosion  suppressant  presents  an  additional  problem.  The  foam 
itself  acts  as  a  secondary  charge  generating  surface  for  the  fuel 
(Ref.  1) .  The  relaxation  time  of  the  charge  in  the  tank  increases 
enormously,  thus  leading  to  substantially  increased  levels  of  accumulated 
charge  in  the  tank,  since  the  charge  is  unable  to  relax  to  earth. 


age 


The  present  study  is  an  attempt  to  apply  numerical  and  computational 
techniques  to  the  problem  in  order  to  examine  the  feasibility  of 
providing  useful  working  guides  for  estimating  electrostatic  potentials 
and  fields  in  fuel  tanks  containing  such  foams. 

THE  MATHEMATICAL  PROBLEM  AND  THE  FINITE  ELEMENT  METHOD 

One  method  of  approach  to  the  problem  is  to  consider  the  fuel  tank 
at  various  levels  of  filling,  postulate  a  charge  distribution  and 
boundary  conditions,  and  solve  Poisson's  equation  for  this  situation. 

This  approach  has  been  tackled  analytically  for  some  very  simple 
geometries  and  charge  distributions  (Refs.  2,3).  Fbr  the  modelling  of 
realistic  situations,  however,  the  analysis  becomes  intractable,  and 
recourse  to  numerical  methods  becomes  essential. 

Various  numerical  techniques  for  approximating  the  solution  of 
electrostatic  field  problems  are  currently  in  use.  These  include  the 
finite-difference  method,  the  charge  simulation  method,  the  boundary 
integral  method,  and  the  finite  element  method.  The  method  chosen  for 
this  study  is  the  finite  element  method.  It  allows  the  modelling  of 
complicated  geometries,  inhomogeneous  charge  distributions,  and 
dielectric  changes  within  the  region  of  interest.  Its  disadvantages 
are  relatively  large  data  input  and  large  computer  storage  requirements. 

The  basic  mathematical  problem  to  be  tackled  is  the  solution 
of  Poisson's  equation 


within  a  region,  subject  to  certain  boundary  conditions.  Here,  4>  is 
the  electrostatic  potential,  p  is  the  space  charge  density,  e  is  the 


absolute  permittivity  (8.854  x  10  12  F/m)  and  C£  is  the  relative 
dielectric  constant  of  the  medium.  To  provide  a  complete  definition 
of  problem,  one  cr  more  of  the  following  boundary  conditions  are  required: 


(i)  4>  =  f(s)  which  fixes  the  potential  at  the  boundary  s,  as  a 

specified  function  f(s)  of  position. 

(e.g.  on  an  earthed  boundary  <p  =  0) 


(ii) 


94> 

3n 


0 


which  forces  equipotentials  to  cross  a  boundary 
normally 


(iii) 


=  h(s)  which  superimposes  a  surface  charge  density 
distribution  on  a  boundary. 


Let  us  now  consider  the  functional 

F=/-j{|V$|2-  <£}ctv  +  /h<f>ds 

v  *  eoer  S 

where  V  is  the  volume  of  the  region  of  interest  contained  by  the 
bounding  surface  S.  If  we  minimise  this  functional,  i.e.  look  for 
admissible  potential  functions  <(>  such  that 

6f  =  0 

it  may  be  shown  that  the  ensuing  Euler- Lagrange  equation  is 


i.e.  the  potential  function  must  satisfy  Poisson's  equation,  subject 
to  appropriate  boundary  conditions. 


Briefly,  the  finite  element  method  utilizes  the  above  results 
by  subdivision  of  the  region  of  interest  into  a  collection  of  elements, 
and  approximating  the  potential  by  a  set  of  piecewise  continuous 
functions  on  these.  Each  element  has  a  fixed  number  of  nodes,  and  the 
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minimisation  is  performed  in  each  element  with  respect  to  the  potential 
values  at  these  nodes. 


For  a  full  three-dimensional  problem,  however,  the  resulting  set 
of  linear  equations  in  the  nodal  potentials  to  be  solved  is  generally 
very  large,  and  necessitates  substantial  computing  resources.  For  many 
problems,  however,  including  the  one  under  consideration,  the  problem 
size  may  be  reduced  in  size  and  complexity  by  using  special  features 
of  the  geometry.  It  will  be  seen  that  for  the  fuel  tank  under  consider¬ 
ation,  the  geometry  may  be  taken  to  be  axisymmetric,  yielding  results 
which  are  sufficiently  accurate  to  make  the  approximation  acceptable. 

To  exploit  axisymmetry,  cylindrical  coordinates  are  adopted,  and 
Poisson's  equation  becomes,  in  a  usual  notation, 


k±  (r  94)  +  1 

r  3r  ,r  3r  1  ? 


+  _ p_ 

Fz*  e  e 

o  r 


where  r,  0,  z  are  cylindrical  coordinates  as  shown: 


T  (r,0,z) 
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For  axisymmetric  problems,  the  geometry  and  charge  distributions 


are  rotationally  symmetric,  the  potential  hence  rotationally 
symmetric,  and  therefore  independent  of  0.  Poisson's  equation  thus 
reduces  to 


o  r 


The  problem  is  thus  basically  two-dimensional: 

I 

Height  from 
Base 

r 

Radial  distance  from  axis  of  symmetry 

The  finite  element  method  consists  of  four  basic  steps.  Firstly  a  grid 
of  numbered  nodal  points  is  established  over  the  region  of  interest, 
including  the  boundary.  At  each  interior  node  the  value  of  potential  is  to 
be  determined.  On  the  boundary  <J>  or  its  normal  derivative,  is  given. 
Secondly  the  nodes  are  interconnected  to  form  a  finite  number  of  subregions 
which  collectively  approximate  the  region  of  interest. 

Thirdly  the  potential  is  approximated  by  a  continuous  function 
over  each  subregion,  continuity  conditions  being  imposed  at  subregion 
boundaries.  Finally  the  unknown  values  of  the  potential  at  the  node 
points  are  calculated  by  means  of  a  variational  principle,  e.g. 
minimisation  of  electrostatic  energy. 

The  subregions  used  in  this  case  are  triangular.  A  first  order 
interpolating  polynomial  (3)  is  used  to  approximate  the  potential 
distribution  over  each  triangle. 

<J>  =  +  a2r  +  a3z  (3) 
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Thus  in  a  triangle  whose  nodes  are  labelled  i,  j,  k  we  have  - 


4>.  = 

a. 

+ 

a„r. 

+  a,z. 

Ti 

1 

2  i 

3  i 

4>.  = 

a. 

+ 

a  r  . 

+  a,z . 

3 

1 

2  3 

3  3 

ii 

-©• 

a. 

+ 

a  r 

+  a„z. 

.1 

2  k 

3  k 

The 

variables 

Ct,  < 

1  2 

<J>.  <j) .  (j>  and  the  coordinates  of  the  noes  i,j,l.  Substituting  back 
i  ]  k 

into  the  interpolating  polynomial  (3)  gives  - 


<j>  =  N.<|>.  +  N  .(J) .  +  N,  <j>, 

11  ]]  kk 


(4) 


where  N. ,  N.,  N,  are  functions  of  r. ,  r.,  r,  ,  z. ,  z.,  z.  ,  r  and  z. 
i  3  k  i  3  k  i  3  k 


8(f>  3d) 

and  may  also  be  readily  obtained. 


The  unknown  values  of  may  now  be  obtained  by  minimising  the 
functional  - 

F  -  f  {r(~)  +  r  (f^)  -  }*dA  -  2  /  $h(s).ds 

A  9r  3z  £o£r  C 

i.e.  6f  =  0 


It  is  not  difficult  to  show  that  this  minimisation  is  equivalent 
to  the  requirement  that  <j>  satisfies  the  axisymmetric  Poisson  equation 
with  appropriate  boundary  conditions. 


.  8d>  8(b 

Substitution  of  the  appropriate  local  expressions  for  <p,  -s^-  and 
into  the  functional,  and  minimisation  by  a  Rayleigh-Ritz  or  equivalent 
technique  yields  a  set  of  simultaneous  equations  in  the  nodal  potentials 


for  each  triangle.  These  elemental  equations  are  then  assembled  to  give 
a  global  set  of  equations  which  may  be  solved  using  standard  techniques. 


3.  DESCRIPTION  OF  COMPUTATIONAL  MODEL 

(i)  Geometry 

The  computational  model  used  is  based  on  the  drawings  of  am  A-10 
fuel  tank  (Ref «  4) .  The  real  tank  configuration  (Figure  1)  is  such  as  to 
permit  an  axisymmetric  approximation.  The  axisymmetric  configuration 
chosen  is  shown  in  Figure  (2) ,  cut  away  to  show  the  inlet  nozzle  and 
the  explosion  suppressant  foam  blocks.  A  dimensional  cross-section 
of  the  tank  is  shown  in  Figure  (3).  (Dimensions  are  in  millimetres). 

(ii)  Boundary  Conditions 

Figure  (4)  illustrates  the  assumed  boundary  conditions.  The  tank 
walls  and  the  inlet  nozzle  are  assumed  to  be  at  earth  potential 
(<fr  =  0) . 

Conforming  to  the  real  situation,  the  foam  is  divided  into  two 
categories,  fixed  and  removable.  Section  4  in  Figure  (4)  corresponds 
to  the  fixed  foam  region.  Sections  1,2,  and  3  correspond  to  removable 
foam  target  sections.  Thus,  possible  target  configurations  are 
Section  1  only.  Sections  1  and  2  together,  or  Sections  1,2,  and  3  together, 
allowing  alteration  of  the  voiding  volume. 

The  dielectric  boundaries  of  the  foam  blocks  are  also  shown  in 
Figure  (4),  as  is  the  axis  of  symmetry.  Since  the  problem  is  axisymmetric, 
the  potential  $  must  satisfy  -  n  on  the  axis  of  symmetry,  i.e.  the 
equipotentials  must  cross  the  axis  normally. 

(iii)  The  Charge  Density 

Accurate  estimations  of  the  volume  and  surface  charge  densities 
occurring  when  fuel  is  pumped  at  speed  into  a  receiving  tank  are  difficult 
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Earthed  Boundaries 


Dielectric  boundaries 
Axis  of  symmetry 

Target  section  1 
Target  section  2 
Target  section  3 
Fixed  foam  region 


Section  of  Tank  Model  showing  Target  Sections 


to  obtain.  Factors  such  as  fuel  conductivity,  rate  of  flow,  turbulence, 
all  contribute  to  the  distribution  of  charge  in  the  fuel  at  any  particular 
time.  The  simplest,  and  not  unrealistic,  distribution  to  assume  is  a 
constant  charge  density  within  the  fuel.  There  is  evidence  (Ref.  5) 
that  jet  fuels  such  as  JP-4  generate  static  electricity  when  passing 
through  poroiis  media.  The  explosion  suppressant  foam  acts  as  a  secondary 
static  charge  generating  surface.  Furthermore,  the  movement  of  charge 
with  foam  present  is  very  slow,  so  a  constant  charge  density  at  any 
filling  level  is  a  reasonable  approximation.  The  dielectric  constant 
of  both  fuel  and  foam  has  been  taken  to  be  2.  A  typical  charge 
distribution  considered  s  shown  in  Figure  (5).  As  we  shall  see,  the 
charge  distributions  studied  at  various  filling  levels  allow  a  wide 
choice  of  postulated  charge  distribution. 


4.  RESULTS 

The  cases  studied  consist  of  charge  distributions,  both  volume  and 
surface,  for  various  filling  levels,  and  for  different  target  sections. 

A  'standard'  constant  volume  charge  distribution  of  10  ** C/m3  in  the 
fuel  was  chosen,  together  with  a  'standard'  surface  charge  distribution 
of  10  C/m  on  the  foam  surface  in  the  voiding  region.  This  is  not 
restricted  since  the  two  charge  distributions  are  treated  separately, 
and  the  principle  of  superposition  allows  simple  scaling  to  a  desired 
charge  density.  The  plots  provided  give  ten  equipotentials  between 
minimum  and  maximum  potentials  for  each  situation.  Also  given  are  plots 
of  electrostatic  field  along  the  axis  of  symmetry  for  each  configuration. 

\ 
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PLOTS  1-17 


The  configuration  for  plots  1-17  consists  of  foam  Section  4  together 

with  foam  Section  1  inserted.  (See  Figure  6). 

Plots  1-5  Charge  density  in  fuel  =  10  4  C/m3 

Charge  density  on  foam  surface  =  0  C/m2 

Filling  levels  at  .1,  .2,  .3,  .4,  .5  metres  above  base. 

Plot  6  Surface  charge  density  of  10  3C/m2  on  Section  1 

upper  surface  only. 

Plots  7-11  Surface  charge  on  Section  1  upper  surface  and  at 

heights  .1,  .2,  .3,  .4,  .5  metres  above  this  surface. 

Plot  12  Vertical  stream  of  fuel  with  charge  density  10  4  C/m3 

impinging  on  the  target  area. 

Plots  13  -  17  Vertical  stream  of  fuel  +  voiding  region  filled  to 
heights  .1,  .2,  .3,  .4,  .5  metres. 
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PLOTS  18  -  30 


The  configuration  for  plots  18  -  30  consists  of  foam  Section  4  together 
with  foam  Sections  1  and  2  inserted.  (See  Figure  7 ) . 


Plots  IS  -  22  Charge  density  in  fuel 


“  10_4  C/m3 


Charge  density  on  foam  surface 


0  C/m 


Filling  levels  at  .1,  .2,  .3,  .4,  .5  metres 


above  base. 


Plot  23 


Surface  charge  density  of  10  3  C/mz  on  Section  2 


upper  surface  only. 


Plots  24  -  26  Surface  charge  on  Section  2  upper  surface  and  at 
heights  *1,  .2,  .3  metres  above  this  surface. 


Plot  27 


Vertical  stream  of  fuel  with  charge  density 


10  4  C/m3  impinging  on  the  target  area. 


Plots  28  ~  30  Vertical  stream  of  fuel  +  voiding  region  filled 


to  heights  .1,  .2,  .3  metres. 
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PLOTS  31  -  39 


The  configuration  for  plots  31-39  consists  of  foam  Section  4  together 
with  foam  Sections  1,  2,  and  3  inserted.  (See  Figure  8). 

Plots  31  -  35  Charge  density  in  fuel  =  10  4  C/m3 

*  ^ 

Charge  density  on  foam  surface  *  0  C/m 

Filling  levels  at  .1,  .2,  .3,  .4,  .5  metres  above  base 

Plot  36  Surface  change  density  of  10  3  C/mz  on  Section  3 

upper  surface  only. 

Plot  37  Surface  charge  on  Section  3  upper  surface  and  at 

height  .1  metre  above  this  surface. 

Plot  38  Vertical  stream  of  fuel  with  charge  density  10  4  C/m3 

impinging  on  the  target  area. 

Plot  39  Voiding  region  filled  with  fuel. 
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APPENDIX  A 


PROGRAM  VALIDATION 

The  finite  element  program  used  in  this  study  was  written  in 
FORTRAN.  To  test  the  accuracy  of  the  program,  it  was  checked  against 
a  problem  for  which  an  analytical  solution  has  been  found.  Asano  (Ref.  6) 
has  provided  an  analytical  solution  for  the  potential  in  a  cylindrical 
earthed  metal  tank  half  filled  with  charged  liquid,  assuming  a  uniform 
charge  density.  This  problem  was  solved  using  the  finite  element 
program  with  a  mesh  of  760  nodes,  and  an  equipotential  map  was  drawn 
(Figure  Al) •  The  parameter  chosen  for  comparison  was  the  field  along 
the  cylinder  axis.  As  can  be  seen  from  Figure  A2»  the  finite  element 
program  produced  results  in  good  agreement  with  those  of  Asano. 
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